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Abstract 

We investigate the behavior of a dilute quasi two-dimensional, harmonically 
confined, weakly interacting Bose gas within the finite-temperature Thomas- 
Fermi approximation. We find that the thermodynamic properties of the 
system are markedly different for repulsive and attractive interactions. Specif- 
ically, in contrast to the repulsive case, there appears to be a phase transition 
when the atoms interact with an attractive pseudo-potential, in the sense that 
there is no self-consistent solution for the normal ground state below a cer- 
tain temperature T* . These numerical findings are supported by analytical 
investigations of the thermodynamics of the system in the complex fugacity 
piane, and within the random-phase approximation. We also show that the 
temperature T* can be interpreted as the limiting temperature below which 
the system cannot be described as a collection of noninteracting haldons. 

PACS: 03.75.Fi, 05.30.Jp, 31.14.Bs 
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I. INTRODUCTION 



Over the last decade or so, there have been a number of theoretical studies of Bose- 
Einstein condensation (BEC) in inhomogeneous low-dimensional Bose systems 
principle, such low-dimensional systems can be created in the laboratory through a suitable 
manipulation of the potential used to trap the ultra-cold atoms. In this paper, we will be 
primarily interested in the situation for which the confinement in the z-direction is much 
steeper than the in-plane confinement, so that the system can be viewed as being (quasi) 
two- dimensionai (2D) in the sense that the motion of the atoms in the ^-direction is "frozen 
out" . 

It is well-known that a homogeneous 2D Bose gas has no long range order (at T 7^ 0), 
and so a BEC transition at any (non-zero) finite temperature cannot occur. However, in 
a 2D harmonically confined system with oscillator frequency uq, it is easily shown that an 
ideal Bose gas has a criticai BEC temperature 

= {G/Tr^y/^N'/'hcoo, (1) 

which, for ^ 1, is much larger than the oscillator gap hujQ. It is interesting to consider 
the limit N 00, coq 0, such that N^^'^ujq = Constant. In this situation, the criticai BEC 
transition temperature Tj°) remains the same. Note that this is not the usuai thermody- 
namic limit, which in 2D would demand that Nujq = Constant, resulting in no BEC. We 
have previously shown that a self-consistent Thomas- Fermi (SCTF) model [0 (i.e., with no 
condensate density), has a solution at ali temperatures provided the mean-field contribution 
of a repulsive zero range interaction is included, no matter how weak. This result contradicts 
a previous claim in the literature |Q that no self-consistent solution of the TP equation in 
a trap is possible below a certain criticai temperature Te. (It is worth pointing out that in 
Ref. [0, the absence of a solution below Te was interpereted as a phase transition in the 
sytem, but not a phase transition to the BEC state.) The absence of a phase transition 
for the repulsive case was also confirmed analytically by examining the branch-points of 
the density in the complex fugacity piane. Note, however, that the SCTF method cannot 
address the important issue of phase fiuctuations in the condensate [^. It was also found 
that the such a system could be described by noninteracting haldons obeying the generalized 
fractional exclusion statistics (FES) |jlO[- This has been recently confirmed by Hansson et 
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For a dilute Bose gas in a three dimensionai trap, where the average distance between 
the atoms is much larger than the scattering length, the Gross-Pitaevskii (GP) approach 
of using a zero-range pseudo-potential with a strength given by AiTh^a/m, (where a is 



the s-wave scattering length, and m is the mass of the atom) is very successful. For a strictly 
2D system however, the strength has a logarithmic energy dependence [0. Recently, Petrov 
et al. P] have deduced the strength of the quasi-2D system by considering free two-particle 
scattering in the piane of the gas, and strong harmonic confinement in the direction (i.e., 
uq -C uJz)- The effective two dimensionai pseudo-potential may be written as a zero-range 
momentum-dependent interaction with a strength given by 



-9 = — - ^ Hnpl,il)]-' , (2) 



m m J2tx 
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where pi2 = (pi — P2)/2 is the relative momentum between two bosons. In the above 
equation, m is the mass of the atom, a the s-wave scattering length, and 1^ = \J^/ 'mojz ■ 
Note that the kinetic energy of relative motion, E = pl2/'m, reduces to twice the single- 
particle energy in the frame where the CM momentum P = (pi + P2)/2 = 0. To simplify 
actual calculations in the many-body problem with this coupling strength, the momentum- 
dependence in Eq.(||) is replaced P by a density dependence by putting = = 2m|/i|, 
where n is the chemical potential. For the bosons in a harmonic oscillator trap, the spatial 
density is inhomogeneous, and the locai chemical potential /i(r) = /i — |mc<jQr^ will be used. 
For the situation when 3> a, the first term on the right-hand side (RHS) of Eq. 
dominates over the second, and the logarithmic term may be neglected. This yields a 
"Constant" interaction strength, 



9 = a/4 , (3) 

m m 

which agrees with our definition for the coupling Constant in an earlier work 0. 

In this paper, we perform finite-temperature SCTF calculations with the above density 
dependent strength factor in the pseudo-potential. Both the repulsive, and the attractive 
cases are considered for realistic values of the physical parameters. For the repulsive case 
(^f > 0), the earlier conclusion that there is no BEC at any temperature remains unaltered. 
For the attractive case {g < 0) on the other hand, we find compelling numerical and an- 
alytical evidence suggesting the existence of a phase transition at a T* which is greater 
than the criticai temperature for BEC given by Eq. (|l]). The outline for the rest of the 
paper is as follows. In Sec. II, we present numerical results of the SCTF calculations for 
both g and g as defined above, as well as some analytical arguments explaining the origin 
of the criticai temperature T*. In Sec. Ili we provide additional support for our numerical 
results by analytically investigating the thermodynamic properties of the system at zero and 
finite-temperatures. Finally, in Sec. IV, we present our concluding remarks. 



II. THOMAS-FERMI CALCULATION 

A. Repulsive Interaction 

We consider bosonic atoms in a two-dimensional harmonic oscillator trap above the 
criticai temperature, if any. In this subsection, a repulsive interaction is used, while in the 
next subsection the interaction is taken to be attractive. In the SCTF method []^, the 
inhomogeneous number density of the interacting bosons is given by 

nir) = I (4) 

J [exp[(é + \mu;lr^ + U{n{r)) - - 1] ' 

where U{n{r)) = ^^^gn(r) is the mean-field potential |]14| generated by the zero-range two- 



body pseudo-potential with the Constant strength g, and is the chemical potential that 
serves to fix the total number of atoms in the trap, 

N = J n(r) dh. (5) 
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To obtain n(r) for a fixed N, Eqs. @ and have to be solved self-consistently. (In our 
numerical work, we have scaled ali lengths by the characteristic in-plane oscillator length 
io = ^Jh/muo, and ali energies by huQ.) The case for ^ > 0, has already been addressed in 
an earlier paper where it was shown that the self-consistent solution for n(r) could be 
obtained for any T > 0, indicating that there is no BEC. We first want to check if the same 
conclusion holds when g is replaced by the more realistic momentum dependent g, as given 
by Eq. (H). As noted before, we replace the pi2— dependence in g by an effective density 
dependence. Under this assumption, the p— integration in can be done analytically. At 
zero-temperature, the result is given by 

^^^n(r) = (/i - ^rnu^r^) . (6) 
m 2 

The RHS of Eq. is precisely the locai chemical potential /i(r). Hence we replace the 
on the RHS of Eq. (0) by {2nhy gn{r) and g itself is determined self-consistently by Eqs. 
and d^). At finite-temperatures, we continue to follow the same procedure, but in this case 
g{r) and n{r) must be found self-consistently from Eqs. (^) and (^). Note that Petrov et al. 
IQ] in one of their footnotes recommend an identical prescription. Following this method, we 
stili find that for a fixed boson number A^, Eqs. and (^) may be solved self-consistently 
for ali T > 0, indicating that there is no phase transition for a repulsive interaction. The 
behavior of the chemical potential /i as a function of T/Tj°^ is shown in Fig. 1. The lowest 
solid curve, for noninteracting bosons, shows the usuai discontinuity in fx at T/Tj:^^ = 1. 
The next higher pair of closely spaced curves show the results for A^ = 10^ Rb^^ atoms 
in a trap with with hujQ = 0.6 nK, Uz = lO^cJo, and a = 5.8 nm. The solid curve shows 
the result for g = 0.05, obtained from the above physical parameters, while the slightly 
lower dashed curve is for the corresponding g given by Eq. (^). We note that the difference 
between the two is negligible for the parameters chosen above, which suggests that FES is 
stili obeyed to a good approximation for relatively small values of g. In the same figure, 
the top three curves are for an artificially boosted interaction, with the continuous curve 
given by g = 0.25. Note that this could in principle be achieved in a variety of ways by 
altering both iz and a. Two limiting cases are given by keeping iz Constant and decreasing 
a by a factor of fi ve, or keeping a fixed and increasing iz by a factor of fi ve. Although both 
procedures yield the same g, we see from Eq. (|]) that g gets affected more when £z alone 
is altered. In Fig. 1, the lowest (dashed) curve is for this case, while the next higher one 
is for the case where only a is decreased by a factor of five. We see that in ali cases with a 
repulsive interaction, the main conclusion that there is no strict BEC continues to hold. 

Even though there is no BEC in the sense that the SCTF solution is obtained without 
paying any special attention to the lowest single-particle quantum state, it is found that 
there is a significant enhancement in the peak of the density distribution as the temperature 
is reduced below Tj:^\ the criticai temperature for the noninteracting case. This is shown 
in Fig. 2 for g = 0.25 (continuous curves) at T/T^^^ = 0.8 and 1.1. When the density- 
dependent version of the coupling g is used, the enhancement is even more pronounced. 
The reason for this behavior is clear from the inset, where the effective g{r), as obtained 
from the self-consistent solution, is shown as a function of r. We note that for the density- 
dependent case, the interaction g{r) is considerably weaker than the Constant g, thus causing 
less suppression of the centrai density. 
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B. Attractive Interaction 



We first take the example of Li^ atoms for which the s-wave scattering length is a = —1.45 
nm. For hujQ = 7 nK (the corresponding oscillator angular frequency uq in the piane is 
271 X 145 rads/s ) and Uz = IO^uq, the couphng Constant g = —0.012. On the other hand, 
the interaction strength could be considerably larger for Cs^'^^ atoms, which have a = —60 
nm. In Fig. 3, we display for = 10^ atoms of Li^ the behavior of the chemical potential n 
as a function of T obtained self-consistently. No solution is found for T < T*, irrespective 
of whether one uses g or g. The highest lying solid curve is for g = —0.012, and difference 
between this and the density-dependent form of g for this weak interaction is too small to 
be shown in the same scale. In the same figure, we also display the effect of boosting the 
attractive interaction strength nearly ten-fold to ^ = —0.1. Not only does the temperature 
T* increases significantly, but it is also possible now to differentiate the curves for the 
cases where the density-dependent form of g is used rather than g. Furthermore, when the 
boosting of the coupling g is done by altering iz rather than a, T* is larger. Nevertheless, the 
behavior of the chemical potential with temperature, for T > T* does not alter significantly 
with the strength of the attractive interaction. Comparing Fig. 3 with the corresponding 
repulsive cases (see Fig. 1), we note that /i showed a much larger variation, but only in the 
range T < Tj°\ which is in any case inaccessible for the attractive interactions. 

Some insight as to why no self-consistent solution for attractive interaction strength g 
may be found by performing the p— integration in Eq. analytically (valid for either sign 
of ^ ). The result is given by 

(yU — -mt^Qr^) = g'n{r) + — ln[l — exp{—2TTh'^n{r)l3/m)] . (7) 

For a negative g, both terms on the RHS are negative, and therefore /i, which is independent 
of r, must be negative at ali temperatures. For this case, if we take the above expression 
literally at T = 0, the density is given by 

= O l~lf.2 (l^l + o"^^0^ ) • (8) 

2'K\g\h 2 

This solution for n{r) is rejected since it keeps increasing with r ! We demand a normalized 
density n{r) that is a monotonically decreasing function of r, which is only found for T > T*. 
By differentiating Eq. with respect to r, we find that for dn/dr < for ali values of r, 
we must have 

[exp(— — n(r)/?) - 1]"^ - |^| > . (9) 

This may be satisfied only for T > T*, which is related to the centrai density at this 
temperature n(0) by the relation 

T* = ^ — -— . (10) 

ln(l + 1/1^1) ^ ^ 

It is interesting that the above equation may also be derived using the the generalized 
exclusion statistics of Haldane [|10|- It has been shown 0JTl| that when the interaction 
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strength of the zero-range pseudo-potential is taken to be a Constant (which we denote by g 
in this paper), the interacting boson system may be mapped on to a noninteracting system 
of haldons, whose statistica! occupancy factor of a single particle state i with energy is 
given by 

Vi = {w. + 9)~\ (11) 

which must be > 0. Note that Wi obeys the relation 

wfXl + w,y-^ = exp[ie,- fx)P] . 

In what follows, we go to the limit coq ^ 0, N ^ oo, such that N^^'^coq = Constant. We then 
identify the centrai density n(0) of the harmonic oscillator potential as the uniform density in 
this limit. The number density n{0) of these haldons at the center of the harmonic oscillator 
potential is given by 

, ^ , ^ + Wq 
n(0) = ^ n -, (12) 



where Wq defìnes the occupancy of the lowest quantum state through Eq. (pH]). We now use 
these relations for the attractive interaction strength g = ~\g\. Note from Eq. ( |Tll) that for 
the occupancy rjo to be positive definite, wq > \g\, and that the occupancy tends to infinity 
when Wq = \g\- Substituting this limiting value of wo in Eq. (0) at T*, we immediately 
obtain Eq. (0). We should regard T* as the limiting temperature below which the system 
cannot be described clS ci coUection of ideal haldons. 

In Fig. 4, the self-consistent spatial density n(r) is shown for = 10^ Li^ atoms in the 
quasi 2D trap, with g = —0.012. For this very weak interaction, the corresponding curve 
for the density-dependent coupling strength g cannot be differentiated on this scale. We 
see that not only is there there a marked enhancement in the centrai density when T is 
decreased from 1.1 T* to T*, but the shape of the density distribution is also different. We 
cannot, of course, get any solution for T < T*. It may be instructive to fìnd a relation 
between T* and analogous to Eq. (|l|), which will of course also depend on the interaction 
strength g. Empirically, we find that this is given by 




= 7(1^1) 77^ àu^o, (13) 



where 7(|^|) = [1 + f (|^|)^''^]- This relation is quite accurate even for small and large g, 
as depicted in Fig. 5. In this figure, we show the numerically calculated points by circles, 
and the continuous curves are derived from Eq. ([T3|). Although the solid curves fitting the 
data have been extended down to A^ = at T = 0, we were not able to numerically access 
stable values of A^ for T* ^ 8, and so Eq. ( [T3| ) is only used to fit the data down to a finite 
temperature T* ^ 8. Also note from Eqs. (|l]) and ([T3|) that T* /TjP'^ is independent of A^, 
and therefore we can draw a "universal" curve of g versus T* /Tjp\ which is shown by the 
inset in Fig. 5. The normal bosonic phase is above the continuous curve shown in this inset, 
but the SCTF method used cannot access the region below. 
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III. ANALYTICAL CONSIDERATIONS 



In this section, we wish to address the question that naturally arises from the numerical 
work presented in Sec. IIB; namely, are our numerical calculations suggesting that there is 
phase transition in the quasi-2D trap for attractive interactions ? Since we are unable to 
answer this question within the framework of the SCTF model, we now turn our attention 
to various analytical approaches which may aid in our assessment of the numerical results 
presented above. 

A. Zero- Temperature Variational Calculation 

Let US consider the 2D analogue of the GP energy functional, namely 



E= (Fr 



(14) 



where g is given by Eq. (0), and may be of either sign. In Eq. (0), the function il){r) 
represents the wave function of a system for which there is macroscopic occupation of the 
lowest quantum state (i.e., ali N bosons are in the ground state): 




^2 



^(0 = i\/-exp(-^) , (15) 
where h is taken to be a variational parameter. We find that the energy is a minimum for 



h=J^{l + N~gY" . (16) 
V mojQ 



For an attractive interaction, g = —\g\, Eq. ([T6| ) implies that the condensate coUapses for 

N\~g\>l. 

This rough estimate is in good agreement with the numerical estimate of A^|^| > 0.94 by 
Adhikari which was carefuUy obtained through a self-consistent solution of the 2D GP 
equations. For a weakly interacting gas like Li^ with g = —0.012, we find that the condensate 
is unstable for iV > 83. 



B. Finite- Temperature Thermodynamics 

This subsection, deals with the study of the behavior of the branch points of the quantity 
u = 27rh'^ (3n{r) / m in the complex fugacity piane Z = exp(/3/i) as a function of the coupling 
Constant g. This is analogous to the analysis that was made by Sutherland while studying 
the thermodynamics of a one-dimensional gas interacting with an inverse square potential 
I3| . The important point to be taken from Sutherland's work is that a branch point of u 



on the positive real axis of Z indicates a phase transition in the system. 
We begin by noting that Eq. (|^) can be rewritten as 
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U 



(17) 



As before, going to the limit — ^ oo, cuq — > with N^l'^uji^ = Constant, Eq. ( P^D reduces to 
the r-independent relation 



_ e(9-'i^)u^ 



(18) 



We note that for ^ = 0, the above equation gives u = — ln(l — Z), implying that u has 
a branch point at Z = 1, i.e. for ^ = 0. This is precisely where BEC takes place for a 
noninteracting gas in a harmonic trap. For g = 1/2, the branch points of u are ai Z = ±2i, 
and for g = 1, the branch point is at Z = — 1. For arbitrary values of g, the branch points 
of u can be obtained by examining the analytical structure of the solutions to the equation 
{x = exp(u)) 



— X 



9'^-Z = 0. 



(19) 



FoUowing Sutherland's prescription ||T5[, we obtain the branch points of u, denoted by Zq, 
for ali ^ > 0: 



(i) ^ > 1 



9Ì9 



~9-l 



(9-1/2) 



(20) 



(ii) 1 > ^ > 



~\ (9-1/2) 
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(21) 



In Fig. 6, the loop and the line along the negative real axis — 1 < Z < in the complex 
fugacity piane shows the branch points of u for positive values of g given by Eqs. ( PU] ) and 
Since there is no branch point on the real fugacity axis for ^ > 0, there is no phase 



transition, in agreement with our numerical calculations. 

The situation for ^ < can be addressed by either solving Eq. 
that the transformation [B!C 



9 



-9 + 1, 



u 



directly, or by noticing 



(22) 



leaves Eq. ([T8|) unchanged. More specifically, under Eq. (p2D, the branch points of u on the 
negative Z— axis between Zq G [—1,0] for 1 < ^ < oo are now shifted to the positive side 
of the axis, namely, Zq G [0, 1] for ali — oo < ^ < 0. The branch points for < ^ < 1 are 
ali mapped back onto Eq. (pi]), thereby leaving the loop structure in the complex Z-plane 
invariant under (^) (see Fig. 6). In other words, for every ^ < we have an associated 
branch point on the positive real Z-axis, suggesting that there is a phase transition in the 
system for attractive interactions. It is important to note, however, that the above analysis 
can offer us no insight into the precise nature of the phase transition; is it BEC ? Is it the 
Kosterlitz-Thouless transition 



Is it the coUapse of the gas 
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C. Random-Phase Approximation 



In this subsection, we will try to address the nature of the phase transition discussed 
above by applying the well-known random-phase approximation (RPA) to investigate 
the stabihty of the trapped 2D Bose gas with attractive interactions. Our approach parallels 
the recent investigation by MeuUer and Baym with the main difference being that we 



confine ourselves to strictly two-dimensions. The basic idea behind their method is to 
identify the collapse of the gas with a divergence in the density-response function x{k,uj; (3), 
which is a measure of the response of the system to some external probe with wave vector 
k and frequency u. As in Ref. |T^, x is calculated within the local-density approximation 



(LDA), which amounts to replacing the response function of the inhomogeneous system by 
that of a system of uniform density n. In order to make contact with the trapped gas, the 
uniform density n is then taken to represent the centrai density of the trap n{r = 0). 

For a confined Bose gas, the collapse is associated with an instability of the lowest energy 
"breathing" mode of the system. In 2D, the breathing mode is given by 



Sp{r) « - l) ^M-r'/il) 



where r is the 2D radiai coordinate, and io = ^h/mujQ is the characteristic size of the cloud. 
The Fourier transform of óp{r) has a maximum at A; = 2/io, and so we look for an instability 
at this wave vector. Owing to the fact the breathing mode has a vanishing frequency, the 
line of collapse is characterized by 

x{k = 2/eo,u; = 0-X)-' = 0, (23) 

where T° denotes the temperature for which the gas undergoes collapse. Note that need 
not coincide with the BEC transition temperature Te of the interacting gas. Since our SCTF 
calculations deal only with the normal phase of the Bose gas, we will restrict our attention to 
the noncondensate contribution to x{k, ^'i T) above, and present a more detailed analysis of 
the RPA in 2D attractive Bose gases elsewhere |^ . Within the RPA, the response function 



for the interacting Bose gas has the structure |21 



X''''\k, u; T) = ^ , (24) 

where Xo the "bare" noncondensate response of the noninteracting system. In 2D, we 
have 

n.. Ti W ^2 /(g - k/2) - fjq + k/2) 
^°(^'"'^)=(2;0^i "^^.-ieig + km-eig-km ' ^'^^ 

In Eq. (^5]), f{q) is the Bose distribution function, and e{k) = h^k"^ /2m is the free particle 
kinetic energy. After integrating out the angular dependence, we have 
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Using the above expression for Xo "we can obtain the hne of collapse by setting k = 2 /io 
and solving for the roots of the denominator in Eq. (^4|). This procedure yields the hne of 
cohapse in terms of /x and T^pa- FoHowing the arguments made in Ref. ||19[, we can then 




relate fi io the total number of normal phase atoms, N, in the harmonically confìned gas 
through the relation 

(27) 

where g^iz) = ^pZ^/p^ is the polylogarithm function. Note that evaluating Eq. ( p^ ) at 
/i = defines the line of condensation for a noninteracting 2D Bose gas. (This result can 
easily be derived by integrating Eq. (^ with ^ = over ali space.) 

The results of our calculations for g = —0.012, —0.1, —0.2 are shown in Fig. 7. For 
comparison, we have also included as dashed lines in the figure the g = —0.1, —0.2 SCTF 
results from Sec. IIB (see also Fig. 5). This figure suggests that the line of collapse derived 
from the RPA has the same functional relationship as the empirical formula given by Eq. (|13D 
in Sec. IIB. In order to confimi this observation, we recali that the relationship used to link 
the uniform gas results to the harmonically confined gas, viz., Eq. (|27|) , is strictly valid only 
for a noninteracting system. As we saw in Sec. IIB, the attractive interactions have the 
effect of "renormalizing" the ideal N vs. T dependence (see Eq. (|l3l)), which will result in 
different curves for the RPA and SCTF data. With this in mind, we rescale Eq. ( P7D by 
1/7(1^1), and comapre the results with the SCTF data. After this procedure, we find that 
the SCTF and RPA curves are almost indistinguishable on the scale of Fig. 7. In this sense, 
we can identify T^pa with the T* of Sec. IIB. This analysis leads us to believe that the 
phase-boundary generated from the SCTF in Sec. IIB is in fact the line of collapse of the 
Bose gas as one approaches the instability from above the BEC transition temperature Te. 

To dose this section, we mention that we have also performed RPA calculations for the 
3D interacting Li'^ Bose gas with attractive interactions (exactly as in Ref. [|1^] but without 
exchange) and found that the 3D gas is more stable than its 2D counterpart; that is, for the 
same temperature and particle number, the 3D system is always in the stable uncoUapsed 
phase relative to the 2D gas. To illustrate this, we have included as an inset to Fig. 7, the 
lines of collapse (solid curve) and condensation (dashed curve) for the 3D system, focusing 
only on the and T regimes that are relevant for a comparison with the corresponding 2D 
system. 



IV. SUMMARY 

One of the first results of this paper was a confirmation of the fact that, at least within 
the SCTF method, there is no strict BEC in a quasi-2D trap for atoms with repulsive 
interactions, even if the more "realistic" moment um-dependent pseudo-potential recently 
derived by Petrov et al. [§, and given by Eq. (H) is used. Indeed, our results indicate 
that the much simpler density-independent coupling Constant g (Eq. (^) is quite adequate 
for numerical work on 2D Bose systems. This observation provides further support to the 
conclusion obtained earlier [0 using a momentum-independent coupling Constant g. 

In the case of attractive interactions, we have found strong numerical and analytical 
evidence for the existence of a phase transition in the system. Since the qualitative features 
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of this phase transition are not strongly influenced by using a more complicateci momentum- 
dependent pseudo-potential, it is again sufficient to use the simpler g in numerical work. 
Our numerical investigations reveal tliat the phase transition is signaled by a "criticai" 
temperature T* > TJ°^ below which there is no self-consistent solution for the normal ground 
state of the system. From the point of view of FES, the T* was interpreted as the limiting 
temperature below which the system cannot be described as a coUection of noninteracting 
haldons. 

In order to further our understanding of the numerical results, we considered several 
analytical approaches. This began in Sec. IIIA with a zero-temperature 2D variational 
calculation of the GP energy functional. Surprisingly, this simple analysis revealed that 
if there is a BEC for attractive interactions in 2D, then the region of stability for the 
condensate (i.e., prior to coUapse) is severely restricted by the condition A^|^| > 1. For Li^, 
this implies an instability for N > 83. The T = calculation was followed in Sec. IIIB by 
an investigation of the branch points of u in the complex fugacity piane. We found that for 
every ^ < 0, there is an associated branch point of u on the positive real Z-axis. This result 
strongly supports our numerical calculations which indicate that there is a phase transition 
taking place in the system for attractive interactions. 

Finally, in Sec. UIC, an RPA-type calculation was used to investigate the line of collapse 
for a 2D Bose gas with attractive interactions. The purpose of this study was to try and 
determine the nature of the phase transition suggested by the analysis in Sec. IIIB. We 
found that the line of instability in the RPA calculation had the same vs. T dependence 
as the SCTF curves in Sec. IIIB. The correspondence between the two formulations is highly 
suggestive that the phase transition is in fact associated with the collapse of the 2D gas as 
one approaches the instability from above Te. We also performed analgous RPA calculations 
for the 3D Bose gas with attractive interactions, and found that for the same and T, the 
3D system is always in a stable, uncoUapsed phase relative to the 2D gas. We will present 
a more detailed comparison of the 2D and 3D RPA calculations elsewhere. 
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FIGURES 



FIG. 1. Chemical potential /x as a function of T/Tc ■ The parameters used to generate this 
figure are consistent with those of the Rb^'' experiments. Sohd hnes are for dashed hnes are for 
g{r) with fixed a and variable (.z, while long-dashed-short-dashed hnes are for g{r) with fixed (.z 
and variable a. Unless stated otherwise, this convention is used in ali subsequent figures. 

FIG. 2. Calculated self-consistent density profile for Constant g = 0.25 and an artificially 
boosted density-dependent interaction g{r). The lowest set of curves correspond to T/Tc^^ = 1.1 
and the highest set to T/Tc^'^ = 0.8. The inset illustrates the self-consistent g{r) evaluated at 
T/xjP^ = 0.8. Notice that g{r) always lies below the corresponding Constant g value. 

FIO. 3. Chemical potential /x as a function of T for attractive interactions. Open circles denote 
the last point for which a self-consistent solution to Eqs. (^) and ^ could be obtained while the 
arrows serve to clarify the location of T* . The parameters used to generate this figure are consistent 
with those of the Li'' experiments. The same line labeling convention as in the repulsive case has 
been used. See text for details. 

FIO. 4. Calculated self-consistent density profile for an attractive interaction g = —0.012. The 
lowest curve is for a temperature T = I.IT* and the highest lying curve is for a temperature 
T = T*. Li^ parameters have been used. 

FIO. 5. Total number of normal phase atoms A'^ as a function of the criticai temperature T*. 
From left to right, each curve is evaluated at g = 0.0, —0.1, —0.2, —1.0. Open circles denote numeri- 
cai data, and the continuous curves are obtained from Eq. (|l^). Inset: Universal "phase-boundary" 
separating the normal bosonic phase from a region which is inaccessible within the SCTF approx- 
imation. 

FIG. 6. The thick continuous curves represent the location of branch points of u in the complex 
fugacity piane as a function of g. Branch points of u lying to the positive real Z-axis indicate a 
phase transition in the system. 

FIG. 7. Total number of normal phase atoms as a function of temperature T. The solid curves 
(RPA), from left to right, are for g = -0.012,-0.1,-0.2 and correspond to N vs. Trpa- The 
dashed curves (SCTF), from left to right, are for g = -0.1,-0.2 and correspond to N vs. T*. 
Inset: Lines of collapse (solid) and condensation (dashed) for the 3D Li^ Bose gas with attractive 
interactions. In this inset, we have zoomed in to show only the N and T ranges relevant for a 
comparison with the corresponding 2D system. See Ref. [19] for details. 
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